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I.  INTRODUCTION 


Consider  the  two  dimensional  vector  ordinary  differential  equation: 

v  =  A  v  +  sgn  (rv1  +  sv2  -  c)  ,  (1) 

2  2 

where  r  +  s  >  0  ,  is  some  non-zero  constant  vector,  and  A  is  a 

constant  matrix  whose  characteristic  polynomial  is 

A  +  2yA  +  (y  +  go  ),  y>0 

Equations  of  this  type  commonly  arise  in  systems  involving  ideal  relays. 

We  refer  to  the  curve  rv^  +  sv^  -  c  =  0  as  the  switching  line  in  these 
applications,  and  to  the  cases  c  ^  0  and  c  =  0  as  biased  and  unbiased, 
respectively.  In  this  paper,  we  investigate  necessary  and  sufficient 
conditions  under  which  (1)  will  have  periodic  (limit  cycle)  solutions. 

We  shall  show  that  existence  of  a  limit  cycle  in  the  unbiased  case  is  a 
necessary  and  almost  sufficient  condition  for  existence  of  a  limit  cycle 
in  the  biased  case,  however,  that  not  all  biasings  (i.e.,  values  of 
c  ^  0)  will  produce  cycles.  We  also  develop  an  extremely  powerful  and 
interesting  geometrical  interpretation  of  the  necessary  and  sufficient 
conditions,  and  use  this  interpretation  to  demonstrate  a  bifurcation 
into  two  cycles  for  the  same  value  of  c  under  certain  circumstances. 

Equations  of  this  type  have  been  considered  by  Fleishman  [2]> 

it 

Davis  and  Fleishman  [l]>  and  Flugge-Lotz  [3],  among  others.  However,  Davis  and 
Fleishman  consider  primarily  the  special  cases  when  only  even  order 

ti 

derivatives  are  present,  and  Flugge-Lotz  does  not  give  a  complete 
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characterization  of  the  limit  cycles.  Also,  neither  author  considers  the 
question  of  multiple  limit  cycles  in  the  general  case. 

For  the  purposes  of  our  investigation,  it  is  convenient  to  introduce 
a  non-singular  linear  transformation  and  time  scaling  (Appendix  A)  that 
reduces  (1)  to: 

x1  =  -28  -  (82  +  l)x2  +  a1  sgn  (x2  -  c)  , 

x2  =  xx  +  a2  sgn  (x2  -  c)  ,  (2) 

2  2 

where  8  >  0  and  al  +  a2  >  ® 

This  now  represents  a  damped  system,  forced  by  a  relay  switching  along 
the  transformed  switching  curve  =  c  .  It  can  be  shown  that  a 

necessary  and  sufficient  condition  to  avoid  the  phenomenon  called 
"chatter”  in  (2),  and  to  insure  that  a  solution  will  exist  at  each  point 
in  the  (x^,x^)  plane,  is  ^  0  .  We  now  make  that  restriction. 

Furthermore,  (2)  involves  only  odd  functions  of  x^  and  x^  ;  hence  we 
can  assume,  without  loss  of  generality,  that  c  ^  0 


2 


II.  THE  PHASE  PLANE 


Our  analysis  of  (2)  begins  with  the  observation  that,  viewed 
separately  in  the  regions  >  c  and  <  c  ,  (2)  reduces  to  a  linear, 

constant  coefficient  system.  Elementary  methods  then  yield  the  solutions 


: k  *  — g  £ 

x^(t)  =  -  c^e  (8  cos  t  +  sin  t)  +  €£0  (cos  t  -  8  sin  t)  -  , 


*  *  -Rr  *  -Rt  (2Ba  +  a  ) 

x  (t)  =  c  e  cos  t  +  c„e  sin  t  +  - r - 

1  (6  +D 


(3a) 


in  the  region  x2^*-^  >  c  *  an<^ 


~  B  t  ”  3 1 

x  (t)  =  -c  e  (6  cos  t  +  sin  t)  +  c0e  (cos  t  -  6  sin  t)  +  a  , 


2 

(3b) 


-St  ,  ,  -St  (26a2ta1) 

x0(t)  =  c  e  cos  t  +  c9e  sin  t  -  - r - 

1  (8^  +  1) 


for  X£(t)  <  c  .  Thus  the  trajectories  in  the  phase  plane  are  simply 
(arcs  of)  decaying  (since  8  >  0)  counterclockwise  spirals,  centered  at 


(- 


(26a  +  a  ) 

- - - — —  ) 

(62  +  1) 


for  the  region  x^  >  c  ,  and  at 


(  a 


2 


(26a  +  a  ) 

- - - —  ) 

2  •  ' 

(6  +  1) 


(4) 


(5) 
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for  the  region  <  c  .  These  centers  will  be  asymtotically  stable 
critical  points  if  and  only  if  they  are  actually  above  and  below 

=  c  ,  respectively,  e.g.,  (4)  describes  a  critical  point  if  and  only  if 

26a2  +  a 

- 5 -  >  c  . 

(ez  +  1) 

Note  that,  if  and  only  if  c  =  0  ,  the  origin  is  also  a  critical  point, 
but  obviously  an  unstable  one* 

A  representative  view  of  the  phase  plane  behavior  can  be  obtained  by 
considering  the  region  x2  >  c  .  There,  in  the  case 

(2Ba2  +  a  )  <  c(B2  +  1)  , 

the  center  for  trajectory  arcs  is  located  in  the  region  <_  c  y  with 

an  x^  co-ordinate  of  -cl^  0  •  Thus,  because  the  spiral  arcs  are 

counterclockwise,  every  arc  emanating  from  a  point  in  x2  >  c  intersects 

the  switching  line  at  an  x^  <  0  in  some  finite  time,  and  exits  the 

region.  Furthermore,  since  the  angular  frequency  is  unity,  then  each 

arc  spends  no  more  time  than  tt  in  the  region,  with  a  time  equal  to  tt 

occurring  only  when  the  center  of  the  arcs  is  located  on  the  switching 

2 

line  (i.e.,  when  (2£a2  +  a^)  =  c(8  4*  1)  ) . 

The  situation  when 

(2Ba2  +  ap  >  c(g2  +  1) 

is  significantly  different  in  that  the  spirals  are  now  centered  at  a 
stable  critical  point  actually  located  in  the  region  x^  >  c  .  Thus, 
some  trajectories  never  exit  the  region,  but  decay  down  to  the  critical 
point,  while  others  intersect  the  switching  line  and  exit  the  region. 
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The  "boundary"  between  the  arcs  that  exit  the  region  and  those  that  decay 
to  the  critical  point  is  the  arc  that  terminates  on  the  switching  line, 


directly  below  the  critical  point,  i.e.,  at 

(  * 


(6) 


This  arc,  which  also  exits  the  region,  is  denoted  (*)  in  Figure  1.  We 


let  T  denote  the  time  represented  by  this  arc,  where  clearly 


7T  <  T  <  2tt  .  It  is  fairly  easily  shown  that  every  other  arc  which  exits 
o 


the  region  x0>  c  represents  less  time  in  the  region  than  T  *  ,  i.e., 

/  o 

* 

T  is  the  maximum  time  an  arc  can  spend  in  x,>  c  and  still  exit  there, 
o  2 

Tq  can  be  determined,  using  (x^  (t),  x^  (t))  to  denote  the 
solution  in  the  region  x^  >  c  ,  by  noting  that  the  condition  described 
by  equation  (6)  and  curve  (*)  of  Figure  1  is  eauivalent  to 


*1*0  ’  '“2 

x  (T  )  =  x  (0)  =  c  . 
z  o  z 


(7) 


Using  (3)  and  simplifying  the  resulting  three  equations,  we  arrive  at: 

* 

-3T  *  * 

e  °  =  cos  T  -8  sin  T  ,  (8) 

o  o 

which  can  easily  be  shown  to  have  one  and  only  one  solution  in  (tt,  2tt)  , 
* 

hence  T^  is  uniquely  defined. 

The  analysis  of  the  region  x^  <  c  is  similar,  and  so  is  omitted. 

The  only  point  of  interest  in  this  case  is  that,  when,  a  critical  point 

2 

is  located  in  this  region  (i.e.,  +  a^)  <  c(8  +1))  a  "maximum 

time",  T  can  be  defined,  similar  to  (7)  by: 


x  (T  )  =  a 
1  o  2 


x  (T  )  =  x  (0)  =  c  . 
z  o  z 


(9) 


5 


S traighf orward  calculation  shows  satisfies  (8),  hence,  since  the 

solution  is  unique, 


T  =  T  *  (10) 

o  o 

(For  the  remainder  of  this  discussion,  we  drop  further  reference  to 
* 

T  ,  replacing  it  by  T  .)  Again  we  can  show  that  no  curve  which 
o  o 

exits  x^  <  c  can  spend  more  than  time  there. 
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Ill,  THE  CONDITIONS  FOR  A  LIMIT  CYCLE 

To  derive  the  necessary  and  sufficient  conditions  for  the  existence 
of  a  limit  cycle  solution  to  (2),  we  begin  with  the  solutions  as  given 
by  (3).  We  let  T^  denote  the  time  spent  by  one  trajectory  arc  in 
<  c  ,  and  T^  the  time  on  an  arc  in  x^  >  c  .  Then,  since  the 
system  (2)  is  autonomous,  and  since  trajectories  must  be  continuous  in 
the  (x^,x^)  plane,  it  follows  that  a  necessary  and  sufficient  condition 
for  a  limit  cycle  of  period  (T^  +  T^)  to  exist  is: 


X1 (Tf)  =  xx  (0) 


X1 (°)  =  \  (t2) 


x2(°)  =  x2^Ti^  =  c 

x  *(0)  =  x  *(T  )  =  c  ^ 


T  ,  T  <  T  . 

1  2  —  o 


(11) 


where  it  is  understood  x^(0)  represents  x^(0  )  >  etc. 

The  condition  given  by  (11)  can  be  reduced  by  using  the  solutions 
given  by  (3)  and  simplifying  the  resulting  expressions  by  eliminating 


c2  and  c2 


This  leads,  in  the  case  where  neither  T^  nor  T2 


equal  tt  ,  to  the  equations: 


cosh  6T-  -  cos  T-  . 

c  - i - 1  +  c 

1  sin  T„  1 


cosh  BT2  -  cos  T2 
sin  T0 


=  0 


sinh  BTj 
'1  sin  T 


-  c. 


*  sinh  BT2 
sin  T„ 


=  -2 


(1-6  )a2  -  6a1 

(62  +  1) 


Tl^To 
T2  —  Tq 


(12) 


(13) 


where 


*  (2Ba?  +  a1) 

c  -  c  =  2  Z  1 

C1  1 


2  ’ 

(6  +  1) 


(14) 
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and  the  value  of  c  for  which  this  cycle  occurs  is  given  by 


c 


=  2  (C1  + 


(15) 


These  equations  express  the  necessary  and  sufficient  conditions  for 

existence  of  a  limit  cycle,  except  in  the  cases  noted.  (Note  these  later 

cases  correspond  respectively  to  c^  =  0  ,  if  and  only  if,  =  tt  and 

* 

c^  =  0  ,  if  and  only  if,  =  tt  .)  In  these  cases,  expressions  similar 

•k 

to  (12)- (14)  can  be  derived.  For  example,  in  the  case  c^  =  0  ,  c^  ^  0  , 
the  condition  for  existence  of  a  limit  cycle  of  period  (T^  +  tt)  is: 


sinh  3tt 

cosh  $7T  +  1 


cosh  gT2  -  cos 
sin  T2 


sinh  0T  (1-0  )ot  -  20a 

- L  +  - - - —  t  <  x 

sin  T2  (20a2  +  a  )  ’  2  -  o  ' 


(16) 


We  now  develop  an  interesting,  and  powerful,  geometric  interpre- 

2  *2 

tation  of  the  necessary  and  sufficient  conditions  when  c^  +  c^  >0 

Consider  (12- (14),  which  represent  the  necessary  and  sufficient  conditions 
* 

when  c^c^  f  0  ,  and  let 


F(t) 


sinh  0t 

sin  t 


G(t) 


cosh  0 1  -  cos  t 

sin  t 


(17) 


Then  (12)  -  (14)  is  equivalent  to  the  vector  equation 
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G(Tp 

-g(t2) 

F(T1) 

* 

"C1 

F(T2) 

1 

1 — 1 

M  _ 

-  — 

2(26a2  +  a^) 
(6 2  +  1) 


(18) 

0 

(1  -  B2)a2  -  Ba1 
(2Ba2  +  ax) 

1 


if  ( 23 +  a  )  ^  0  .  But  this  means  that  the  three  vectors  in  (18)  are 

3 

linearly  dependent,  i.e.,  lie  in  the  same  plane  in  R  .  Clearly  this 
plane  cannot  coincide  with  the  z  =  1  plane,  yet  since  each  vector  in 
(18)  terminates  on  the  z  =  1  plane,  linear  dependence  then  requires  that 
the  end  points  of  the  three  vectors  lie  in  a  line  in  that  plane.  Thus, 
(12)-(14)  can  be  satisfied  if  and  only  if  the  vectors 


G(T  ) 

-G(T2)- 

0 

9 

,  and 

(1  -  62)  a  -  6a- 

F(T1) 

f(t2) 

Z  1 

-(26a2  +  a1) 

(19) 


lie  in  the  same  straight  line.  But  by  drawing  the  parametric  curves: 

?!  =  { (G(x) ,  FCc))  0  <  x  <  T  }  ,  (20) 

1  —  —  o 

C2  =  { (-G  (x)  ,F(t))  0  <  x  <  Tq} 

plus  the  point 

(1  -  B2)  a  -  6a 

P  =  (0, - - - -  )  ,  (26a  +  a  )  4  0  ,  (21) 

0  (26a2  +  ap 

'fc 

we  have  that  (provided  ^  0  ’  anc^  (26a2  +  a^)  4  0)  ,  any  line 
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through  P  and  both  curves  generates  a  limit  cycle,  and,  conversely, 
o 

for  any  limit  cycle,  there  must  be  a  straight  line  through  the  point  and 
both  curves. 

Straightforward  but  laborious  calculation,  plus  liberal  application 

of  hyperbolic  identities,  shows  that  both  F(t)  and  G(*r)  are  positive 

and  monotonically  increasing  on  (0 , tt )  ,  and  negative  and  convex  down  on 

(tt,  2tt)  .  Furthermore  G(t)  has  its  maximum  to  the  right  of 

(i.e.  G(t)  increases  monotonically  on  (tTjT^))  ,  while  F(t)  has  its 

maximum  at  a  T  <  T  .  Also,  and  are  convex  up  in  the  upper 

mo  ±  z 

half  plane,  intersecting  at  (0,3  )  ,  and  convex  down  in  the  lower  half 
plane,  being  asymptotic  to  the  lines 


and 


sinh  3^  0 

U  cosh  3tt  +  1  w 


(22) 


.sinh  Btt  _ 

cosh  3tt  +  1  W 


(23) 


respectively,  as  t  tt  (w-*-±°°)  .  Thus  and  must  have 
the  approximate  shapes  shown  in  Figure  2.  Figure  3  shows  one  situation 
which  would  lead  to  a  limit  cycle.  Observe  from  (12)  that  the  ratio  of 
the  horizontal  intersections  on  and  >  which  we  denote: 

* 


G(TX) 

g(t2) 


(24) 


satisfies  (from  (14)  and  (15)); 


c 


1 

2 


(C1  + 


(23a2  + 

(3  2  +1) 


(1  +  P) 


(1  ~  p) 


(25) 
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Thus,  since  p  is  just  the  ratio  of  horizontal  co-ordinates,  we  could 
easily  compute,  from  the  geometrical  configuration,  the  value  of  c  for 
the  cycle  shown  in  Figure  3.  (Assuming  values  for  B^  ,  and 

are  known.) 

Equations  (24)— (25)  will  be  especially  useful  in  discussing  the 

range  of  values  of  c  for  which  cycles  can  occur,  and,  since  it  will 

become  important  later,  note  that  on  any  interval  not  including  p=  1, 

c  as  given  by  (25)  is  a  monotonically  increasing  [decreasing]  function 

of  p  for  (2go<2  +  a^)>  0  [(230^  +  a  )<  0] 

In  the  discussion  so  far,  this  geometrical  interpretation  has  been 

limited  to  the  cases  c^c^  ^  0  >  and  ^Ba^  +  a^)  ^  0  .  However,  the 

2  *2 

interpretation  can  easily  be  extended  to  all  cases  where  c^  +  c^  >0 
by  extending  the  interpretation  of  a  line  intersecting  the  point  and  both 
curves  to  allow  intersection  in  the  limit  at  infinity.  Thus,  for  example, 
the  case  c  =  0  ,  c^  ^  0  and  (230^  +  a^)  ^  0  is  treated  (refer  to 
(24))  as  the  limiting  case 

G(T  )  +  ±  »  . 

Now,  referring  to  Figures  2-3,  it  is  obvious  “intersection”  in  this 
limiting  case  occurs  if  and  only  if  the  straight  line  is  through  Pq  , 
parallel  to  the  asymptote  of  £  ,  intersecting  ^  •  (Figure  4).  That 

is,  there  must  be  a  T^  such  that 


(1  -  B2)  a  -  6a 

F(T2>  +  ■  (26^  +-a-2r 

-  G(T2) 


sinh  Btt _ 

cosh  Btt  +  1 


(26) 


But  referring  to  the  definitions  of  F(t)  and  G(t)  it  can  be  seen  (26) 
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is  identical  to  (16).  A  similar  process  can  be  followed  when 
c1*  =  0  ,  c^  ^  0  ,  and  (23a2  +  a ±  0  . 

Lastly,  interpreting  (19)-(21)  where  (23a^  +  a^)  =0  as  the 
limiting  case  when  Pq  approaches  the  point  at  infinity  along  the  verti¬ 
cal  axis,  we  see  that  an  intersection"  can  occur  if  and  only  if  the 
points  on  and  ^  li e  on  the  same  vertical  line,  or  equivalently 

G(T^)  =  -  G(T2)  •  However,  since  by  (14)- (15)  this  case  corresponds  to 

* 

cl=  c±  =  c  ,  this  interpretation  immediately  satisfies  (12),  and  yields 
from  (13)  a  cycle  occuring  for  every  value  of  c  satisfying: 


c 


sinh  3T^ 
sin 


sinh  3T2 
sin  T2 


(27) 


where  and  are  solutions  of  G(T^)  =  ~G(T^)  ,  and  the  absolute 

value  of  the  quantity  in  brackets  represents  the  vertical  distance 
between  the  two  points  of  intersection. 

Thus,  including  the  limiting  cases,  we  see  that  for  all 

2 

cl  +  ci  >  0  ,  the  geometrical  interpretation  is  equivalent  to  the 
solution  of  (12)- (14),  and  hence  it  expresses  the  necessary  and  sufficient 
condition  for  existence  of  a  limit  cycle. 


12 


IV.  EXISTENCE  OF  LIMIT  CYCLES 


In  this  section  we  apply  the  geometric  construction  developed  above 
to  derive  necessary  and  sufficient  conditions  under  which  a  limit  cycle(s) 
will  exist.  Our  main  conclusion  will  be  to  determine  existence  conditions 
in  terms  of  3  >  and  •  Specifically,  we  show  a  bound,  depen- 

dent  solely  on  3  ,  such  that  no  cycles  exist  for  any  values  of  c  if 
the  ratio  [a^/a^]  exceeds  this  bound,  and,  for  all  values  less  than  this 
bound,  there  is  a  c  >  0  ,  depending  on  3  ,  a.  ,  and  such  that 

a  cycle  will  exist  if  and  only  if  c  <  c 

J  J  ii—  max 

We  begin  by  noting  the  following  result,  alluded  to  in  the  introduc¬ 
tion  to  this  paper: 


Lemma  1:  A  limit  cycle  will  exist  for  values  of  c  >  0  if  and 
only  if  a  cycle  of  the  same  period  exists  for  -  c 


The  proof  of  this,  which  follows  straightforwardly  from  the  oddness 
of  all  functions  in  (2),  is  omitted  here. 

The  following  theorem  gives  the  condition  under  which  no  cycles  of 
any  period  are  possible,  irrespective  of  the  biasing. 


Theorem  1: 


There  is  an  n  (3)  such  that  (2)  will  have  no 
max 


limit  cycle  solutions  if  either 


(a)  c*2  =  0  ,  or 


(b) 


and 


>  n  (8) 

max 


Proof:  Let  M(B)  =  -  max 


(tt>  2tt) 


sinhBt 


sin  t 


>  0, 


(28) 
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and 


n 

max 


(8) 


1  -  8  -  28M 

(6  +  M) 


>  -28 


(29) 


(Note  -  M  is  the  maximum  vertical  coordinate  on  both  £  and  5  in 
the  lower  half  plane.)  It  is  then  easily  shown  that  conditions  (a)  and 
(b)  are  equivalent  to 


-  M  < 


(1  -  62)a2  -  30^ 
(28a2  +  ap 


However,  this  implies  that  Pq  (as  given  by  (21))  lies  above  the  line 
joining  the  maxima  of  and  ^  in  the  lower  half  plane,  but  at  or 

below  the  juncture  of  and  ^  in  the  upper  half  plane  (Figure  5). 

The  convexity  of  and  5^  and  the  quandrants  in  which  the  curves 

are  located  immediately  imply  no  straight  line  through  Pq  can  intersect 
both  ^  and  ^  >  hence  no  cycle  can  exist.  | 

Note  that,  by  defining 

a2  4  0  ,  (30) 

Pq  can  be  equivalently  defined 


PQ(B,n) 


(1  -  82)  -  8ti 

(23  +  n) 


(31) 


Then,  if  the  situation  when  =  0  is  considered  as  the  limiting  case 

n  "*■  ±  00  ,  the  previous  theorem  implies  cycles  can  exist  only  if 

—  oo  <  n  <  n  (3)  •  We  now  complete  the  result  by  showing  that  if 
—  max 

”  00  <  n  £  T)  (3)  »  a  limit  cycle  will  exist  for  at  least  one  value  of  c 
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Theorem  2:  If  n  =  n  (B)  ,  a  cycle  will  exist  if  and  only  if 
-  max  J 

c  =  0  . 

Proof:  If  n  =  n  (6)  ,  then 

-  max 

Pq  =  (0,-M)  , 

that  is  P  lies  on  the  line  joining  the  maxima  C-  and  £  in 
o  12 

the  lower  half  plane,  (Figure  5).  But  the  convexity  of  and 

imply  this  is  also  the  only  line  through  P  intersecting  both 
2  o 

curves.  Thus  one  and  only  one  cycle  will  occur  in  this  case.  From 
Lemma  1,  c  =  0  .  | 


The  above  result  is  the  reason  why,  in  the  introduction,  we  commented 

that  existence  of  a  limit  cycle  in  the  unbiased  case  was  a  necessary  and 

almost  sufficient  to  give  cycles  for  some  biasings,  for,  as  the  following 

shows,  only  when  n  (6)  =  n  can  a  cycle  exist  only  for  c  =  0 

max 

Theorem  3:  If  -  00  <  n  <  n  (6)  ,  then  there  is  a  c  >  0  , 

-  max  max 

such  that  a  limit  cycle  will  exist  if  and  only  if 


Proof :  By  Lemma  1  we  restrict  ourselves  to  c  >_  0  .  The  actual 

proof  is  divided  into  three  subcases;  -26  <  n  <  n  (6)  , 

max 

n  =  -  26  ,  and  -  00  <  n  <  -26 

Case  I:  The  case  -26  <  n  <  n  (6)  >  corresponds  (with  our  re- 
-  max 

striction  a ^  0)  to  (260^  +  ot^)  >  0  ,  and  so,  by  (25),  c  >  0 

corresponds  to  -1  <_  p  <  1  .  Also,  from  (29)  -  (31),  we  can  show 

that  in  this  case  P  (n>B)  is  located  in  the  lower  half  plane, 

o 

below  the  line  joining  the  maxima  there  of  anc^  *  However, 
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TT 

with  P  so  located,  we  can  define  0  ,  0  <  0  <  —  ,  as  the 

o  o  o  2 

largest  clockwise  angle  from  the  horizontal  for  which  a  straight  line 

through  P  intersects  in  the  third  quandrant.  Then  clearly 

o  i 

any  line  through  Pq  with  slope,  ,  will  intersect  both  an^ 

C 2  if  and  only  if 

|  m  |  <_  tan  0^  .  (Figure  6) 

However,  since  p  must  vary  continuously  with  m^  (along  lines 

intersecting  both  curves),  and  since  p  can  approach  1  only  if  the 

intersecting  line  approaches  the  vertical,  then  p  must  be  bounded 

away  from  1  ,  and  there  exists  a  p  such  that  intersection  will 

J  max 

occur  for  c  >  0  if  and  only  if  -  1  <  p  <  p  <  1  .  It  follows 

—  max 

from  (25)  that  a  cycle  will  exist  for  c  >_  0  if  and  only  if 


(2ga?  +  a  ) 

0  <  c  <  c  »  — - — 

-  "  (g2  +  x) 


(1  +  p 
(1  -  P 


) 

max 

max^ 


> 


since  the  right  hand  side  of  (25)  is  strictly  increasing  for 

(2Ba„  +  a, )  >0  and  -1  <  p  <  1  .  Clearly  c  depends  on  3  ,  a. 

2  1  —  max  l 

and  <*2 

Case  II:  When  -23  =  n  ,  (23ot2  +  a^)  =  0  ,  and  as  we  have  discussed 
before  in  this  case  limit  cycles  exist  if  and  only  if  c  satisfies 
(27),  where  in  (27),  the  quantity: 

sinh  3T  sinh  3T2 


represents  the  distance  between  £  and  t,  on  any  vertical  line 
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intersecting  both  (since  Pq Cn )  is  imagined  passed  to  a  limit  at 

infinity  along  the  vertical  axis).  Since  the  distance  must  clearly 

have  a  non-zero  minimum,  and  since  (27)  involves  its  reciprocal, 

it  follows  straightforwardly  that  there  is  a  c  such  that  c  , 

max 

as  defined  by  (27)  satisfies 

0  <  c  <  c 
—  —  max 


Case  III:  When  -  °°  <  n  <  -23  ,  the  proof  is  essentially  similar 

to  case  I,  except  P  (h,3)  is  located  in  the  upper  half-plane, 

o 

above  the  juncture  of  and  ^  •  | 


Theorem  3  completes  the  description  of  the  necessary  and  sufficient 
conditions  under  which  a  limit  cycle (s)  will  exist.  Note  that  although 
this  theorem  predicts  at  least  one  cycle  for  every  value  of  c  satisfy¬ 
ing  0  <  c  <  c  ,  it  does  not  address  the  possibility  of  more  than 
—  1  1  —  max 

one  cycle  occuring  at  a  single  value  of  c  .  We  shall  pursue  this  in 

the  next  section.  However,  we  first  investigate  in  some  more  depth  the 

properties  of  n  (3) 
max 

Elementary  methods  show  that,  for  a  given  3  >  M(3)  as  defined  by 

(29)  is  given  by 


where  T  (3) 
m 


M(3)  =  -  ‘ 

is  the  solution,  in 


sinh  3T  (3) 

_ m 

sin  T  (3) 
m 

(tt,2tt)  ,  of  the  equation 


tan  t  =  “  tanh  3t 

p 

Since  the  right  hand  side  of  this  equation  can  be  shown  to  be  monotoni- 
cally  decreasing  in  3  (for  fixed  t  )  ,  it  follows  that  Tm(6)  is 
monotonically  decreasing  in  3  >  and 
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lim  T  (3)  =  tt 

n  V  m 


Further  direct  calculation  will  show 


,m/oN  T  (6)  cosh  3T  (3) 
dM(3)  _  _m _ m 

d3  ~  sin  T  (3) 


>  0 


m 


since  tt  <  T  .  Thus  M(B)  is  a  monotonically  increasing,  and  clearly 
m 

M(B)  00  as  $  ^  00 

With  M(B)  thus  fully  characterized,  it  is  easily  shown  that 

ri  (B)  »  as  given  by  (29),  is  a  monotonically  decreasing  function  of 
max 

B  ,  asymptotic  to  the  curve  r|  =  -2B  as  B  00  .  Figure  7  shows  the 

general  dependence  of  r>max  >  and  it  is  fairly  obvious  that 

n  (3)  =  -23  ,  3  >  1  . 

max 
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V.  BIFURCATION  OF  THE  LIMIT  CYCLES 


The  previous  section  showed  that  for  -  00  <  H  Ji  max(3)  ,  there  was 

a  c  (n,3)  such  that  at  least  one  limit  cycle  would  exist  for  every 
max  J 

c  satisfying  0  <  |c  <  c  .  In  this  section  we  complete  the 

description  of  the  limit  cycles  by  showing  there  are  values  of  c  that 
allow  more  than  one  cycle,  how  these  values  are  determined,  and  how  the 
transition  from  single  to  multiple  cycles  occurs. 

Consider  the  situation  represented  by  figure  8.  When  the  straight 
line  intersects  two  curves  and  ^  >  the  horizontal  co-ordinates  of 

the  intersections  and  slopes  there  of  these  curves  are  denoted 
W1  *  mi  anC*  w2  »  ra2  respectively,  and  a  is  the  angle  between  the 
line  and  the  horizontal.  Then,  using  elementary  calculus,  it  can  be 
shown  that: 


{  —  >  =  —  { 
da  w  ^  w  ^ 


m 


+  1  } 


(m2-m1) 


(m  -m  )  (m  -m  ) 
la  z  a 


However,  referring  to  (24),  this  equation  becomes 


it,  .  2  .  tyV 

da  p  ma  ^  (m  -m  )  (m  -m  ) 

la  z  a 


(33) 


where  m^  is  the  slope  of  the  intersecting  line.  This  equation  allows 
us  to  determine,  relatively  easily,  the  points  where  changes  sign. 

But  from  (25),  c  is  a  monotonically  increasing  function  of  p  for 
(23  +  n)  <  0  ,  and  a  monotonically  decreasing  function  for  (23  +  n)  <  0 
Thus  every  time  changes  sign  (for  fixed  values  of  3  ,n  )  ,  a  bifur¬ 

cation  occurs,  i.e.  two  cycles  exist  for  the  same  values  of  c 

With  this  formulation,  and  our  main  intent  being  to  characterize 
the  algebraic  sign  of  as  given  by  (33),  we  proceed.  As  will  emerge 
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more  clearly  later  in  the  discussion,  there  are  several  special  subcases, 

of  increasing  complexity,  that  must  be  considered.  We  shall  not  do  this 

exhaustively,  but  essentially  look  in  full  detail  only  at  one  of  these 

cases,  then  comment  on  the  others.  In  essence,  the  complex  cases  arise 

when  P  is  located  "near”  (0,-M)  ,  so  that  a  single  line  through  P 
o  o 

can  intersect  (or  twice  in  the  same  quadrant. 

Consider  the  case  where  Pq  is  located  in  the  lower  half  plane, 

and  at  or  below  the  intersection  of  the  tangent  to  at  T  with  the 

1  o 

verticle  axis.  (Figure  9).  Since  we  can  show  m. (T  )  =  -1  ,  this  condi- 

1  o 

tion  can  be  reduced,  using  the  definition  of  T  ,  to 


where 
tan  0 


0  <  20  +  n  £  - 

F(t)  was  defined  in  (17).  Note 
>  1  . 


(1  +  B2) 

2 (F (Tq)  "  B) 

also  this  immediately  implies 


o 

In  this  situation,  due  to  the  symmetry  and  convexity  of  and 

,  no  line  thru  P  can  intersect  either  curve  more  than  once  in  any 
2  o 

quadrant.  Furthermore,  as  noted  in  the  previous  section,  since 

(20  +  n)  >  0  it  follows  cycles  with  c  >_  0  occur  only  for 

-1  £  p  <  1  ,  which  equates  to  lines  thru  Pq  whose  slopes  satisfy 

-  tan  0  <  m  <  0  , 

O  —  OF  — 

where  0  was  defined  in  Theorem  3  and  Figure  6  . 
o 

The  easiest  way  to  arrive  at  the  information  needed  to  determine 

is  by  describing  the  variation  of  m  and  m^  as  a  increases  from 

-0  to  zero.  In  this  case,  observe  that  0  is  determined  by  the  line 
o  o 

thru  P  to  (G(T  ),  F(T  ))  ,  the  end  point  of  .  Thus  as  a 

o  o  o  1 

increases  to  zero,  m^  increases  monotonically  (due  to  the  convexity  of 
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from  -1  to  a  positive  value  less  than 

sinh  Btt  .. 

cosh  3tt  +  1 

Since,  as  noted  above,  tan  >_  1  ,  then  as  a  increases  to  zero, 

m  =  tan  a  must  increase  monotonically  from  -tan  to  zero.  Note 

that  since  a  line  through  P  as  located  here  cannot  intersect 

o  *1 

tangentially  (except  perhaps  at  the  endpoint) ,  it  follows 

m  <  m_  for  -0  <  a  <  0 

a  1  o 

Las t ly ,  since  m  <  -1  at  -0  ,  then  the  line  through  P  at  -0 

o  —  o  o  o 

must  intersect  £ ^  second  quadrant.  (Since  here  asymtotic 

to  a  curve  whose  slope  is  greater  than  -1).  Thus  as  a  increases,  m^ 

will  initially  monotonically  decrease  from  a  value  in  (-1,  0)  to 

sinh  3 tt 
cosh  3tt  +  1  * 

at  which  time  the  intersection  point  transitions  to  the  fourth  quadrant 

(and  also  p  passes  through  zero,  changing  sign),  then  m^  increases 

monotonically  to  a  value  at  a  =  0  which  is  the  negative  of  m^  at 

a  =  0  .  It  thus  follows  m^  <  m^  <  0  in  second  quadrant,  and 

m^  <  <  0  in  the  fourth  quadrant  for  a  <_  0 

If  m^  ,  m^  and  m^  are  now  plotted  based  on  the  above  analysis, 

a  picture  like  Figure  10  will  emerge.  Observe  that  and  must 

intersect  exactly  once  for  -0  <  a  <  0  .  Now  if  we  denote  the  inter- 

J  o 

section  of  m,  and  m^  as  -a,  ,  and  that  of  m,  and  m  as  -a  , 

12b  1  a  t 

where  necessarily  -a,  <  -a  ,  we  can  easily  construct  Table  1  to  give 

b  t 

the  algebraic  signs  necessary  to  interpret  (33). 
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TABLE  I 


-0  <  a  <  -  a, 

o  b 


-a,  <  a  <  -a 

D  t 


-a  <  a  <  0 


ml  ■ 

m 

a 

+ 

+ 

+ 

m2  ' 

m 

a 

+ 

+ 

- 

m2 

+ 

- 

- 

P 

+ 

+ 

- 

d£_ 

da 

+ 

- 

- 

From  this  table  it  is  clear  that  multiple  cycles  (for  the  same  value 

of  c  )  can  exist  in  some  cases.  Furthermore  this  multiplicity  can  be 

of  order  exactly  two,  since  only  changes  sign  once.  (Again  recall 

c  is  a  monotonic  function  of  p  .)  Lastly,  the  multiple  cycles  can 

occur  only  in  a  band  of  values  of  c  whose  lowest  bound  is  strictly 

positive.  In  fact,  due  to  the  sign  of  ^  on  -0  <  a  <  -a,  ,  the  value 

da  o  b 

of  c  at  which  the  multiple  cycles  first  occur  is  precisely  the  one 

found  from  the  value  of  p  for  the  line  with  slope  -0  .  It  should  be 

o 

noted  that  since  the  phase  plane  is  two  dimensional,  when  double  cycles 

occur  at  least  one  will  be  unstable.  This  behavior,  existence  of  a 

single  cycle  for  0  <  c  <  c,  ,  then  two  cycles  for  c,  <  c  <  c  ,  seems 

b  b  —  max 

best  described  by  the  statement  that  a  bifurcation  occurs  at  c^ 

The  situation  when  (26  +  ri)  £  0  can  be  easily  constructed  in  a 
manner  to  the  above,  and  yields  essentially  the  same  result  -  a  unique 
cycle  will  exist  for  0  <  I  cl  <c,  where  c,  <  c  ,  there  for 
c^  <_  |  c |  <_  c  exactly  two  cycles  occur. 

The  other  possible  situations  which  arise  become  progressively  more 
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complicated  from  a  bookkeeping  point  of  view.  The  reason  for  this  lies 

in  the  geometry,  with  the  convex  nature  of  the  curves  involved.  For, 

since  is  convex,  the  curve  formed  by  plus  the  extension  of  the 

tangent  line  at  the  endpoint  (at  T  )  to  the  intersection  with  the 

o 

verticle  axis  is  also  convex.  Thus,  where 

.  ik+jh.  <  26  +  n  , 

2 (F (T  )  -3) 
o 

P  lies  "outside"  this  extended  curve,  and  hence  for  at  least  some 
o 

angles  o  ,  a  line  through  Pq  with  slope  m^  will  intersect  twice. 

Thus  m^  cannot  be  described  by  a  single-valued  function  of  a  ,  and 

the  analysis  becomes  more  complex.  However,  as  long  as  P^  is  located 

below  the  line  joining  the  endpoints  of  an d  ^  t'ie  third  and 

fourth  quadrants,  the  same  basic  behavior  as  was  found  in  the  simple 

cases  above  can  be  shown  to  exist. 

However,  when  P  lies  on  the  line  joining  the  endpoints  of  C-, 
o  1 

and  C  ,  i.e.  when 

1  2 

o  <  26  +  1  *  -  pVt  Z  6  .  (34) 

o 

a  qualitative  change  in  the  behavior  occurs.  For  at  this  point,  as 

Figure  11  shows,  there  are  two  solutions  (points  (l)-(3)  and  (2)-(4))  for 

which  p  =  -1  ,  and  hence  c  =  0  .  Thus,  in  addition  to  the  bifurcation 

at  c^  >  0  ,  there  is  now  a  second  bifurcation  at  c  =  0  .  As  Pq 

continues  to  move  toward  (0,-M)  ,  eventually  a  condition  is  reached 

where  multiple  cycles  exist  at  all  values  of  c 

Thus,  in  conclusions,  we  have  shown  that,  in  all  cases,  there  is  a 

c  such  that  two  cycles  exist  for  c,  <  I  cl  <  c  ,  and  in  certain  cases 
b  b  —  1  1  —  max 
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there  is  also  a 


0  < 


1 

Cb 


-  °b 


such  that  two  cycles  can  exist  for 
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VI .  SUMMARY 


In  this  paper  we  have  considered  limit  cycle  solutions  to 
x1  =  -23x1  -  (32+1)x2  +  sgn  (x2-c) 

x2  =  X^^  +  a2  sgn  (x2~c) 

We  have  derived  the  necessary  and  sufficient  conditions  for  limit  cycle 
solutions  to  exist,  and  have  presented  the  conditions  in  terms  of  an 
equivalent,  easily  visualized,  geometric  condition*  With  the  use  of  this 
geometric  interpretation,  we  then  showed  that,  with  one  exception, 
existence  of  a  limit  cycle  for  c  =  0  is  both  necessary  and  sufficient 
for  cycles  to  exist  for  some  range  of  c  4  0  .  Lastly,  we  have  shown 

that  multiple  cycles,  for  some  values  of  c  ,  are  possible  is  essentially 
all  cases* 
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Appendix  1 


We  consider  the  system 


v  =  A  v  +  sgn (rv^  +  sv^  -  c) 


(A-l) 


where  u  is  a  constant  vector,  r  +  s  >0,  and  A  is  a  real, 

— o  = 

2 

constant  2x2  matrix  whose  characteristic  polynomial  is  X  +  2yX  + 
2  2 

(y  4-  to  ) ,  where  y,  a)  >  0.  First  observe  that  the  non-singular 
transformation 


S  ■  B  v, 


reduces  (A-l)  to 


-1 


w  =  (BAB  )w  +  B  u  sgn 


(a)W2  -  c) 


(A-2) 


where,  because  of  similarily,  BAB  has  the  same  characteristic 
polynomial  as  A,  i.e. 


x2  + 


2y\  +  (y  + 


o>2) 


Since  to  >  0,  both  off-diagonal  terms  in  BAB’*'  must  be  non-zero, 
and  so  BAB'*'  will  have  the  form, 
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BAB 


-1 


11 


21 


12 


22 


where  d2^d  12  *  °*  Thus  4  0,  and  hence  the  transformation 


X  =  c  w  , 


22 


— 


is  non-singular.  This  transformation  reduces  (A-2)  to 


±  =  (CBAB1C1)v 


+  C  B  _u^  sgn  (uy2  -  c) 


(A-3) 


where 


C  (BAB1)  C_1 


r 


(d  2 1  +  d22) 


^dlld22  d21d12') 


Again  by  similarity,  this  matrix  has  the  characteristic  polynomial 

2  2  2 
A  +  2y\  +  (y  +  u  )  . 


By  computing  the  characteristic  polynomial  directly,  we  see 


-2y  =  (d 


11 


d22) 


2 

Y 


+  03 


dlld22  d21d12 
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Thus , 


2  2 

yx  =  -  2Xyx  -  (y  +  <*>  )y2  +  ^  sgn  (<oy2  -  c) 


(A-4) 


y2  =  yl 


+  u  sgn  (toy  -  c) 


where 


u  =  C  B  u 


But,  now  the  time  scaling  oat  =  t  reduces  (A-4)  to 


y1  =  ~2(^)y1  -  ((^-)2  +  1)  (coy2)  +  — •  sgn  ((toy2)  -  c) 


(wy2)  =  yx 


+  u2  sgn  ((toy2)  -  c) 


which,  with  the  transformation 


x.,  = 


yl  •  x2  '  "y2  '  S  ‘  1 


a.  =  —  ,  a„  =  u 
1  to  2  2 


becomes  (2) . 
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